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ABSTRACT 

Five planets are known to orbit the star 55 Cancri. The recently-discovered 
planet / at 0.78 AU (Fischer et a/. 2008) is located at the inner edge of a 
previously-identified stable zone that separates the three close-in planets from 
planet d at 5.9 AU. Here we map the stability of the orbital space between 
planets / and d using a suite of n-body integrations that include an additional, 
yet-to-be-discovered planet g with a radial velocity amphtude of 5 ms~^ (planet 
mass = 0.5-1.2 Saturn masses). We find a large stable zone extending from 0.9 
to 3.8 AU at eccentricities below 0.4. For each system we quantify the proba- 
bility of detecting planets 6 — / on their current orbits given perturbations from 
hypothetical planet g, in order to further constrain the mass and orbit of an addi- 
tional planet. We find that large perturbations are associated with specific mean 
motion resonances (MMRs) with planets / and d. We show that two MMRs, 
3f:lg (the 1:3 MMR between planets g and /) and 4g:ld cannot contain a planet 
g. The 2f:lg MMR is unlikely to contain a planet more massive than ~ 20 M^. 
The 3g:ld and 5g:2d MMRs could contain a resonant planet but the resonant 
location is strongly confined. The 3f:2g, 2g:ld and 3g:2d MMRs exert a stabi- 
lizing infiuence and could contain a resonant planet. Furthermore, we show that 
the stable zone may in fact contain 2-3 additional planets, if they are ~ 50 M^ 
each. Finally, we show that any planets exterior to planet d must reside beyond 
10 AU. 
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1. Introduction 

In a remarkable study, Fischer et al. (2008) have measured the orbits of five planets 
orbiting the the star 55 Cancri, the most planets of any exoplanet system to date. The 
system contains two strongly-interacting, near-resonant giant planets at 0.115 and 0.24 AU 
(Butler et al. 1997; Marcy et al. 2002), a 'hot Neptune' at 0.038 AU (McArthur et al. 2004), a 
Jupiter analog at 5.9 AU (Marcy et al. 2002) and a newly-discovered sub-Saturn-mass planet 
at 0.78 AU (Fischer et al. 2008). Table 1 lists Fischer et al. 's self-consistent dynamical fit of 
the orbits of the five known planets in 55 Cancri. 

The fast-paced nature of exoplanet discoveries can lead to interesting interactions be- 
tween theory and observation. Prior to the discovery of planet 55 Cancri /, several groups 
had mapped out the region between planets c and d to determine the most likely location of 
additional planets. Most studies used massless test particles to probe the stable zone (Barnes 
& Raymond 2004 - hereafter BR04; Jones, Underwood & Sleep 2005; Rivera & Haghighipour 
2007). Test particles are good proxies for small. Earth-sized planets because they simply 
react to the ambient gravitational field. However, they are not good substitutes for fully- 
interacting, real planets. Thus, Raymond & Barnes (2005; hereafter RB05) mapped out this 
zone using Saturn-mass test planets. The stable zone from BR04 and RB05 extended from 
0.7 AU to 3.2-3.4 AU, a region that includes the star's habitable zone (Raymond, Barnes & 
Kaib 2006). The planet 55 Cnc / was discovered by Fischer et al. at the inner edge of that 
stable zone. 

The "Packed Planetary Systems" (PPS) hypothesis asserts that if a zone exists in which 
massive planets are dynamically stable, then that zone is likely to contain a massive planet 
(BR04, RB05, Raymond et a/. 2006; Barnes, Godziewski & Raymond 2008). Although the 
idea behind the PPS hypothesis is not new (see, for instance, Laskar 1996), the large number 
of planetary systems being discovered around other stars allows PPS to be tested directly. 
Indeed, the ~ 1.4 Saturn mass planet HD 74156 d recently discovered by Bean et al. (2008) 
was located in the stable zone mapped out in BR04 and RB05, and with the approximate 
mass predicted by RB05 (Barnes et al. 2008). In addition, most of the first-discovered plane- 
tary systems are now known to be packed (Barnes et al. 2008), as well as ~85% of the known 
two- planet systems (Barnes & Greenberg 2007). The fact that 55 Cancri / hes within the 
stable zone identified in previous work (BR04; RB05) also supports PPS, especially since 
planets e through c are packed, i.e. no additional planets could exist between them. Several 
other planet predictions have been made and remain to be confirmed or refuted (see Barnes 
et al. 2008) - the most concrete outstanding prediction is for the system HD 38529 (see 
RB05). 

Mean motion resonances (MMRs) are of great interest because they constrain theories 
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of planet formation. Models of convergent migration in gaseous protoplanetary disks predict 
that planets should almost always be found in low-order MMRs and with low-amplitude 
resonant libration (Snellgrove et a/. 2001; Lee & Pcalc 2002). This may even have been the 
case for the giant planets in our Solar System (Morbidelli et a/. 2007). On the other hand, 
planet-planet scattering can produce pairs of resonant planets in ~ 5% of unstable systems, 
but with large-amplitude libration and often in higher-order MMRs (Raymond et a/. 2008). 
Thus, understanding the frequency and character of MMRs in planetary systems is central 
to planet formation theory. 

In the context of PPS, 55 Cancri is an important system as it contains many planets, 
but still appears to have a gap large enough to support more planets. Therefore, PPS makes 
a clear prediction that another planet must exist between known planets / and d. In this 
paper we add massive hypothetical planets to the system identified by Fischer et al. (2008) to 
determine which physical and orbital properties could still permit a stable planetary system. 
We focus our search on the "new" stable zone between planets / and d. We also show that 
certain dynamically stable configurations are unlikely to contain a planet because the large 
eccentricity oscillations induced in the known planets significantly reduce the probability 
of Fischer et al. having detected the known planets on their identified orbits, to within the 
observational errors. The orbital regions that perturb the known planets most strongly 
correlate with specific dynamical resonances, such that we can put meaningful constraints 
on the masses of planets in those resonances. Finally, we also use test particle simulations to 
map out the region of stability for additional planets beyond planet d, in the distant reaches 
of the planetary system. 



2. Methods 

Our analysis consists of four parts; the methods used for each are described in this 
section. First, we map the stable zone between planets / and d using massive test planets - 
note that we use the term "test planets" to refer to massive, fully-interacting planets. Our 
numerical methods are described in § 2.1.1. Second, we use massless test particles to map 
the stability of orbits exterior to planet d, as described in § 2.1.2. Third, we use the same 
technique to map several mean motion resonances in the stable region. A simple overview 
of resonant theory is presented in §2.2. Finally, we use a quantity called the FTD - defined 
in § 2.3 - to evaluate the probability of detecting stable test planets. 
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2.1. Numerical Methods 

2.1.1. Massive test planets 

We performed 2622 6-planet integrations of the 55 Cancri planetary system which in- 
clude an additional hypothetical planet g located between known planets / and d. In each 
case, the known planets began on orbits from Table 1 including randomly-assigned mutual 
inclinations of less than 1 degree. Planet g was placed from 0.85 to 5.0 AU in increments 
of 0.03 AU and with eccentricity between 0.0 and 0.6 in increments of 0.033. The mass of 
planet g was chosen to induce a reflex velocity of 5 ms^^ in the 0.92 Mq host star (Valenti 
& Fischer 2005): its mass was varied continuously from ~ 50 inside 1 AU to 120 at 
5 AU. Orbital angles of the planets g were chosen at random. The system was integrated for 
10 Myr using the symplectic integrator Mercury (Chambers 1999), based on the Wisdom- 
Holman mapping (Wisdom & Holman 1991) We used a 0.1 day timestep and all simulations 
conserved energy to better than 1 part in 10^. Integrations were stopped when they either 
reached 10 Myr or if a close encounter occurred between any two planets such that their Hill 
radii overlapped. 

Although 10 Myr is much less than the typical ages of extrasolar planetary systems 
(~ Gyr), for a survey of this magnitude it is impractical to simulate each case for Gyrs. 
Previous N-body integrations of extrasolar planets have shown that 10^ orbits is sufficient 
to identify ~ 99% of unstable configurations (Barnes & Quinn 2004). Moreover, N-body 
models of stability boundaries are consistent with alternative methods, such as the Mean 
Exponential Growth of Nearby Orbits (MEGNO; Cincotta & Simo 2000) or Fast Lyapunov 
Indicators (Froschle et al. 1998; Sandor et al. 2007). For example, 1 Myr N-body integrations 
of the 2:1 resonant pair in HD 82943 (Barnes & Quinn 2004) identified a stabihty boundary 
that is very close to that of a MEGNO calculation (Gozdziewski & Maciejewski 2001). 
More recently, Barnes & Greenberg (2006a), using 1 Myr N-body integrations, derived a 
quantitative relationship between the Hill and Lagrange stability boundaries for the non- 
resonant planets in HD 12661 that is nearly identical to a MEGNO study (Sidlichovsky & 
Gerlach 2008). Therefore, for both resonant and non-resonant cases, 10^ year integrations 
provide a realistic measurement of stability boundaries. 

In Section 4, we performed several thousand additional integrations but with hypothet- 
ical planet g in or near specific mean motion resonances (MMRs) with planet / or d. In each 
case we aligned planet (7's longitude of pericenter vj and time of perihelion with either planet 
f or d unless otherwise noted. Small mutual inclinations (< 1 deg) between the two planets 
were included, with random nodal angles. Each set of simulations focused on a given MMR 
and included test planets of fixed mass with a range of orbital parameters designed to cover 
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the MMR. The number of simulations ranged from 30 (4g:ld) to >1100 (2g:ld) simulations 
per set. Planet g^s mass was constant in each set of simulations but varied by a factor of 2-3 
between sets from the maximum value (RV = 5 ms~^) down to 10-40 M®. We performed 
2-3 sets for each MMR. 

Our results are clearly sensitive to the assumed "true" orbits and masses of planets 
b — f. For this work we have adopted Fischer et al. 's (2008) self-consistent dynamical fit, 
but the observational uncertainties remain large. However, the locations of the MMRs in 
question scale simply with the semimajor axis of planet d or /. The strength of these MMRs 
depends on the mass and eccentricity of planets d or f (e.g., Murray & Dcrmott 1999). The 
eccentricity of planet d is relatively well-known, while that of planet / is weakly constrained. 
Thus, the system parameters that could affect our results are e/, and Mf. Since we 
assumed a small value of e/, any increase would affect the strength of the 3f:2g, 2f:lg, 
and 3f:lg MMRs. If Mf and Ma increase due to a determination of the system's observed 
inclination, then all the resonances we studied will increase in strength. This will tend to 
destabilize planets and also increase the size of chaotic zones. Thus, our results are likely to 
be "lower limits" in terms of the strength of resonances. Despite these potential issues, our 
simulations provide a realistic picture of the (in) stability of each MMR. 

2.1.2. Massless Test Particles 

To give a more complete view of the planetary system, we also tested the stability of 
planets exterior to planet d (5.9 AU). We used massless test particles for these simulations 
because of their smaller computational expense. Test particles were spaced by 0.01 AU from 
6 to 30 AU (2401 total particles), and were given zero eccentricity, zero inclination orbits. 
All five known planets were included with orbits from Table 1, including randomly assigned 
inclinations of less than 1 degree. As in previous runs, we used the Mercury hybrid integrator 
(Chambers 1999) with a 0.1 day timestep and integrated the system for 10 Myr. 

2.2. Theory of Mean Motion Resonances (MMRs) 

For mean motion resonance p-|-q : p, the resonant arguments 9i (also called "resonant 
angles" ) are of the form 

Oi,2 = {p + q)Xi - pX2 - q^i,2 (1) 

where A are mean longitudes, w are longitudes of pericenter, and subscripts 1 and 2 refer to 
the inner and outer planet, respectively (e.g., Murray & Dermott 1999). Resonant arguments 



-6- 



effectively measure the angle between the two planets at the conjunction point - if any 
argument librates rather than circulates, then the planets arc in mean motion resonance. 
In fact, the bulk of resonant configurations are characterized by only one librating resonant 
argument (Michtchenko et a/. 2008). In general, libration occurs around equilibrium angles 
of zero or 180° but any angle can serve as the equilibrium. Different resonances have different 
quantities of resonant arguments, involving various permutations of the final terms in Eq. 
1. For example, the 2:1 MMR (q=l, p=l) has two resonant arguments, and the 3:1 MMR 
(q=2, p=l) has three arguments: 

9i — 3Ai — A2 — ^wi, O2 — 3Ai — A2 — 2^2, and 0^ — 3Ai — A2 — {wi + W2)- (2) 

In Section 4, we focus on the possibility of a hypothetical planet g existing in several 
MMRs in the stable zone between planets / and d. We measure the behavior of planets in 
and near resonance using the appropriate resonant arguments, as well as the relative apsidal 
orientation, i.e., Wg — w^j- 

2.3. The FTD value ("Fraction of Time on Detected orbits") 

We have developed a simple quantity to constrain the location of hypothetical planet g 
beyond a simple stability criterion. To do this, we consider the observational constraints on 
the orbits of known planets h - f (1-sigma error bars from Fischer et al. (2008) are listed in 
Table 1). A stable test planet can induce large oscillations in the eccentricities of the observed 
planets. Systems undergoing large eccentricity oscillations can be stable indefinitely as long 
as their orbits remain sufficiently separated (Marchal & Bozis 1982; Gladman 1993; Barnes & 
Greenberg 2006a, 2007). However, systems with large eccentricity oscillations are less likely 
to be observed in a specific eccentricity range, especially with all planets having relatively 
small eccentricities, as is the case for 55 Cancri. The probability that a hypothetical planet 
g can exist on a given orbit is related to the fraction of time that known planets b — f 
are on their current orbits, to within the observational error bars. We call this quantity 
the FTD ( "Fraction of Time on Detected orbits" ) . If the FTD is small, then it is unlikely 
for planet g to exist on that orbit, because perturbations from planet g have decreased the 
probability of the a/reac??/-mao?e-detection of planets b — f. However, if the FTD is close to 1 
then planet g does not significantly affect the hkelihood of detecting the other planets and 
therefore hypothetical planet could exist on the given orbit. We have calibrated the FTD to 
have a value of unity for the known five-planet system (with no planet g). To perform this 
calibration, we artificially increased the observational error of planet c from 0.008 to 0.013. 
This was necessary simply because the evolution of the five known planets causes planet c's 
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eccentricity to oscillate with an amplitude that is larger than its observational uncertainty, 
such that the FTD of the 5-planet system is ~ 0.65. Thus, we calibrate by artificially 
increasing the uncertainty to roughly match the oscillation amplitude. As the region of 
interest lies between planets / and d, low FTD values are virtually always due to increases 
in the eccentricities of planets / or d. The small change we made to the error of planet 
c does not affect our results, and different methods for calibrating the FTD yield similar 
values. The FTD value therefore represents a quantity that measures the perturbations of 
a hypothetical planet g on the detectability of observed planets b — f, normalized to the 
amplitude of the self-induced perturbations of planets b — f. 

To summarize, regions of high FTD (white in upcoming figures) represent orbits of 
planet g which are consistent with current observations of the system. Regions of low FTD 
(blue or black) represent orbits which significantly decrease the probability of detecting 
planets 6 — / on their observed orbits. Thus, we do not expect an additional planet to exist 
in regions with low FTD. Our confidence in this assertion scales with the FTD value itself 
(see color bar in upcoming figures). We a low FTD value to be below 50%, although this 
choice is arbitrary and much of the dynamical structure of the stable region is revealed at 
FTD values above 0.5. Note that all regions that have an FTD value are dynamically stable 
for our 10 Myr integration. 



3. The stable zone between planets / and d 

Figure [1] shows the stable zone between planets / and d: 984 of the 2622 simulations 
were stable (37.5%). Hatched areas indicate unstable regions, white and grey/blue indicate 
stable zones. The inner edge of the stable zone is defined by orbits that approach within a 
critical distance of planet / (the dashed line denotes orbits that cross those of planets / or 
d). The outer regions of the stable zone are carved by resonances with the ~ 4 Jupiter-mass 
planet d. Virtually no stable regions exist exterior to the 2:1 mean motion resonance (MMR) 
with planet d at 3.7 AU, except for the 3:2 MMR at ~ 4.5 AU (not all test planets at 4.5 
AU in Fig [1] are in resonance because angles were chosen randomly) . Note that the outer 
boundary of the stable zone is more distant than the one mapped in RB04 and BR05 - this 
is due to a decrease in the best-fit eccentricity of planet d, reducing the strength of its secular 
and resonant perturbations. For a given semimajor axis and eccentricity of test planet g, the 
bluescale of Fig [T] represents the FTD, i.e. the probability of detecting known planets b — f 
on their current orbits (see color bar). The dark observationally unlikely areas do not fall at 
random, but are associated with specific dynamical structures within the stable zone. The 
wide, dark band from 1.3-2 AU with e ~ 0.2 — 0.4 are orbits for which secular perturbations 
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from planet g increase the eccentricity of planet / above 0.2. The wide dark dip from 2-2.4 
AU at smaller eccentricities is associated with a secular resonance between planets / and 
g which also increases the eccentricity of planet / above its observational limit. All other 
observationally unlikely (i.e., low FTD, dark) regions are caused by MMRs with planets / or 
d, although some are not clearly resolved in Fig. [1] because the resonance is narrow. There 
is clearly room in between planets / and c for an additional planet; in § 5 we explore the 
possibility that multiple companions might lie in this zone. 



4. Mean motion resonances (MMRs) 

We performed extensive additional simulations to test the stability of parameter space 
in the vicinity of eight resonances in the stable zone - 2g:3f (the 2:3 MMR between planets g 
and /), lg:2f, lg:3f, 4g:ld, 3g:ld, 5g:2d, 2g:ld and 3g:2d. The location of these resonances is 
shown in Fig. [T]and listed in Table 2. Based on our results, we divide the eight MMRs into 
three categories: stable, unstable, and neutral resonances. A stable MMR effectively stabi- 
lizes a given region against secular perturbations (i.e., long-term gravitational perturbations 
far from resonance; see e.g. Murray & Dermott 1999). For example, as seen in Fig. [1], there 
are locations associated with the 3g:2d MMR at ~ 4.5 AU that, although they cross planet 
d's orbit, are stable for long timescales. Conversely, an unstable resonance destabilizes a 
region that would be stable under just secular perturbations. For example, the region at 2.8- 
3.0 AU is well-shielded from secular perturbations, but the 3g:ld MMR at 2.88 AU causes a 
large swath of nearby orbits to be unstable. A neutral resonance is one where a region would 
be stable under secular perturbations, and remains stable with the resonance. Although the 
stability of test planets is not strongly affected by these MMRs, FTD values can be strongly 
affected, which in turn affect the likelihood of detecting a planet in a neutral resonance. 

We see general similarities between different resonances. In many cases there exists a 
small region that can undergo resonant libration - that region is usually confined in a^, Cg, 
and Mg (the mass of planet g) space. Planets in this region undergo regular eccentricity 
oscillations such that their FTD values are usually quite high, i.e. a planet can exist in 
that zone. Just outside a resonant region there often exists a chaotic zone in which planets 
may undergo temporary capture into the resonance. These zones are characterized by large 
irregular eccentricity variations that can eventually lead to close encounters and dynamical 
instability. The instability timescale is shorter for smaller Mg such that these chaotic zones 
are more populated for large Mg. However, given the relatively short 10 Myr duration of our 
integrations, we suspect that these chaotic zones would be cleared out in the system lifetime. 
We also found that stable zones with apsidal libration often exist close to the resonance. 
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4.1. Stable Resonances - 3f:2g, 2g:ld, and 3g:2d 

4.1.1. The 3f:2g MMR 

The 3f:2g MMR is located from 1.02-1.04 AU. Figure [2] shows the outcome of 136 
simulations with planet g in the resonant region, formatted as in Fig. [H Two stable peaks 
extend above the collision line with planet /, at 1.024 and 1.034-1.039 AU. To avoid a close 
encounter and maintain dynamical stability, these planets must be in the 3:2 MMR. Indeed, 
the resonance provides a protection mechanism to maintain stability despite crossing orbits. 
The resonant dynamics prevents close encounters from happening by phasing orbital angles 
in various ways (see section 3 of Marzari et al. 2006) - this is also the case for the 2g:ld and 
3gL2d MMRs. As expected, we find that all planets on the two peaks above the collision 
line undergo resonant libration of 9i = 3Xg — 2Xf — Wg about 180°. In the peak at 1.034 AU, 
resonant orbits extend down to zero eccentricity. However, the resonance associated with 
the peak at 1.024 AU extends down to eg ~ 0.05. Below that limit and for the rest of the 
nearby, low-eccentricity stable zone, test planets are not in resonance with planet / 

Figure [3] shows the evolution of a simulation above the collision line with planet /. 
Libration of 61 about 180° is apparent. In contrast, tOg — Wf and 62 are preferentially 
found near 0° but they do occasionally circulate. If all three angles were librating then the 
system would be in apsidal corotation resonance; Michtchenko & Beauge 2003; Ferraz-Mello, 
Michtchenko & Beauge 2003). The eccentricities of planets g and / oscillate out of phase 
with amplitudes of ~ 0.3. Note that ej therefore exceeds the limits of its observational 
uncertainty, since its nominal current value is ~ with an uncertainty of 0.2. Thus, this 
simulation has a low FTD value of 0.335. 

FTD values for test planets above the collision line are smaller for larger values of Mg. 
However, more than half of resonant configurations have very high FTD values. Therefore, 
a planet as massive as 54 M^ could reside in the 3f:2g MMR, but only at low eccentricity 
(e.<0.2). 

4.1.2. The 2g: Id MMR 

The 2:1 MMR with planet d (2g:ld) is a wide, stable resonance located from 3.6-3.85 
AU, and in some cases extending above the collision line with planet d. Figure H] shows the 
outcomes of our integrations near the resonance. There is a peak of stability from 3.6-3.9 
AU, and a sharp cliff of instability for ag > 3.9 AU. The height of the peak depends on 
Mg-. the stable region extends to higher e for more massive planets. The majority of the 
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stable region in Fig. H] participates in the 2g:ld MMR, i.e. at least one resonant argument 
librates. However, the behavior of different resonant arguments varies with Mg Figure [5] 
shows the stable zone from Fig. H] color-coded by which angle is librating {6i = 2\d — A^cCg 
and 6*2 = 2Arf — Ag — Wd)- The libration of 6i is widespread and covers a large area. In 
contrast, 62 librates only in cases with Mg = 100 M^, at the center of the resonance, right on 
the collision line with planet d. In cases where 6*2 librates, 9i and Wg — Wd also hbrate in a 
configuration known as an apsidal corotation resonance. For lower Mg, the apsidal corotation 
resonance is apparent only in a few cases for Mg = 50 M^. It is interesting that the small 
island of 62 libration for Mg = 100 has very high FTD values, while surrounding areas, 
while still in the resonance, have far lower FTD values (Fig. H]). These high FTD areas are 
shifted to slightly higher eg for Mg = 50 and are in fact unstable for Mg = 20 M®. If 
a planet g exists in the 2g:ld MMR, then it must be localized in both mass and orbital 
parameter space. For large Mg, the planet could be either right on the collision line with 
planet d at ag ~ 3.73 AU and eg ~ 0.5, or in the surrounding region of high FTD that 
extends from 2.6-2.85 AU with eg from 0.1-0.4. The lower-FTD belt that separates these 
two regions has FTD ~ 0.7, so we cannot firmly exclude planets from that region. For 
smaller Mg, only the second region is available, although it reaches slightly higher eg. 



4.1.3. The 3g:2d MMR 

The 3g:2d MMR is the most dramatic example of a stabilizing resonance. The entire 
resonant region is unstable to secular perturbations (See Fig. [1]). Nonetheless, Figure [6] 
shows that there does exist a contiguous stable region here. Moreover, more than half of 
the resonant region has orbits that cross that of planet d. We find that all orbits across the 
collision line with planet d exhibit regular libration of the resonant angle 61 = 3Xd — 2Xg — Wg 
about 0°, although none undergo apsidal libration. For the majority of cases below the 
collision line there is a preferential alignment of 61, theta2, and Wg — xud, but circulation 
does occur. The situation is similar for the three different values of Mg, although a larger 
fraction of systems exhibited stable resonant libration for lower Mg. 

FTD values above the colhsion line are 0.5-0.8 for Mg = 113 M®, 0.8-1 for Mg = 50 M®, 
and 1 for Mg = 20 M^. This suggests that the 3g:2d MMR is unlikely to contain a planet 
more massive than ~ 50 M^ above the collision line. However, just below the collision 
line FTD values are large for all masses so we cannot constrain Mg beyond the stability 
boundaries. 

It is interesting that low-eccentricity test planets are unstable in this region. This 
appears to be due to short-term dynamical forcing from planet d, as the low-e^ region does 
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not participate in the 3g:2d MMR. Planet (i's Hill sphere is very large, ~0.65 AU, such that 
any body exterior to 4.88 AU will cross planet (i's orbit unless a favorable alignment (i.e., a 
resonance) prevents this. For a test planet starting at 4.5 AU, an eccentricity greater than 
0.07 will bring the planet into the orbit-crossing region. Secular forcing from planet d is very 
strong in the region of the 3g:2d MMR, so any planet not participating in the resonance 
will be quickly destabilized. For low-e^ orbits near, but not in, the 3g:2d MMR, encounters 
between planets g and d can occur in less than two orbital periods of planet d. 

4.2. Unstable Resonances — 3g:ld and 4g:ld 

4.2.1. The 3g: Id MMR 

The 3g:ld MMR is not truly an unstable resonance, although Figure [7] shows that a 
large region of parameter space centered on the resonance (at ~ 2.88 AU)0 is destabilized. 
However, a small range of test planets does show evidence of long-term stable libration of one 
of the three resonant arguments for the 3:1 MMR (see Eq. 2). This region is located at 
= 2.86-2.89 AU and Cg < 0.06 (i.e., Cg < e^). In these cases only one argument, ^3, librates, 
whereas 61, 62, and Wg — all circulate. The eccentricities of planets g and d oscillate 
regularly within narrow ranges such that the FTD value of these resonant cases is low. In 
other words, a configuration with planet in 3:1 resonance with planet d is observationally 
allowed, although the resonant region is narrow and restricted to very low eccentricities. 

Figure [H] shows the evolution of two simulations, one in stable resonant libration and the 
other undergoing chaotic evolution including a time spent in resonance. In the stable case, 
the apses of planets d and g are circulating but 6*3 librates consistently with an amplitude of 
60°. In contrast, the chaotic (and ultimately unstable) case undergoes resonant libration of 
9i for 1.5 Myr, during which eg remained confined in a relatively narrow band and Wg — Wd 
librated about anti-alignment (see below). Once the resonance was broken, eg ranged from 
close to zero to above 0.5. At 3.2 Myr, planets g and d underwent a close encounter and the 
integration was stopped. 

There exists a small "island" near the resonance at = 2.85-2.88 AU with eg = 0.15—0.2 
which is stable for long timescales. This island is small but apparent for all three test planet 
masses and in all cases the island has high FTD values, i.e., test planets in this region do 
not strongly perturb the orbits of planets b — f. In this island, the longitudes of pericenter 
of planets d and g librate with low amplitude and eccentricities of both planets also oscillate 



The location of the resonance is shifted shghtly from its nominal value of 2.83 AU by secular effects. 
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with relatively low amplitudes. Thus, this island of low-amplitude apsidal libration has very 
high FTD values. There is another region in Fig. [7] which exhibits low-amplitude apsidal 
libration, with a > 2.88 AU and e ~ 0.06 (note that = 0.063). This region is not distinct 
from surrounding orbits in terms of the FTD value; nonetheless it is strongly localized. It is 
interesting that this libration is so strong on one side of the resonance (i.e., at orbital period 
ratios with planet d of less than 3:1) and nonexistent on the other side of the resonance. 

Test planets near the resonant region (a^ = 2.86 — 2.89 AU, Cg < 0.06) or apsidal- 
libration island (a^ = 2.85 — 2.88 AU, Cg = 0.15 — 0.2) may undergo temporary capture into 
the 3g:ld resonance, i.e. temporary libration of one or more resonant arguments. However, in 
these cases the evolution of the system is typically chaotic such that resonant libration does 
not last for long times. The majority of these cases are unstable on the 10 Myr integration 
period, especially for smaller test planet masses Mg. For larger Mg, stable cases have small 
FTD values and so are observationally unlikely. In addition, we expect such cases to be 
unstable on longer timescales given the chaotic evolution of the system. 

FTD values at large Cg are a function of Mg (see Fig. [7]), as a more massive eccentric 
planet will impart larger perturbations on the other planets in the system. Note that these 
regions do not undergo resonant or apsidal libration. 

We reran the same cases with the apses of planets g and d anti-aligned rather than 
aligned; Figure summarizes the outcome. For anti-aligned apses we see the same instability 
of planets in the resonant region, but no island of apsidal libration was apparent. There also 
existed a few cases undergoing stable resonant libration of 63 in the same region as the 
aligned case (a^ = 2.86-2.89 AU), but only for initial eg = 0. The only other test planets 
that underwent resonant libration were for Mg = 90 at higher eccentricities. As before, 
these cases evolve chaotically and have high FTD values. Such orbits are unstable for smaller 
Mg and likely unstable on longer timescales for Mg = 90 M®. 

The stability limits far from resonance differ between the aligned and anti-aligned sim- 
ulations. In particular, the edges of the resonance occur at lower eccentricities for the anti- 
aligned case (at Cg = 0.3 — 0.35 rather than 0.45-0.5). This appears to be due to stronger 
secular forcing for the cases which are initially anti-aligned. In other words, anti-aligned test 
planets start the simulations in a phase of eccentricity growth and aligned planets start in 
a phase of eccentricity decline. Thus, the long-term median eccentricity of planet g in an 
anti-aligned configuration with planet d is significantly larger than the eccentricity of planet 
g starting in an aligned configuration. Higher eccentricities lead to closer encounters with 
other planets, which is the key factor in determining the stability of a planetary system 
(e.g., Marchal & Bozis 1982; Gladman 1993; Barnes & Greenberg 2006a, 2007). Therefore, 
for a given starting eccentricity, a planet in an anti-aligned configuration will have a higher 
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average eccentricity than for an aligned configuration - this higher eccentricity will bring the 
anti-aligned case closer to instability. So, although the stability limit for aligned and anti- 
aligned cases has the same time-averaged eccentricity, this limit occurs for smaller starting 
eccentricities for the anti-aligned configuration. It is therefore important to note that the 
initial eccentricity is not necessarily a good measure of the typical eccentricity during an 
integration, especially when comparing systematically different orbital angles. 

4.2.2. The 4g:ld MMR 

The 4:1 MMR with planet d is strongly dependent on Mg (see Figure fTU]) . For both 
Mg = 80 M^ and 40 M^, the outskirts of the resonance at high-FTD values show the same 
structure. However, the heart of the resonance, at 2.35-2.36 AU, is populated with lower- 
FTD planets for Mg = 80 Mq and is empty for Mg = 40 M®. Planets in this region undergo 
chaotic and temporary capture into resonant libration. However, the resonance never persists 
for more than a few Myr. For Mg = 40 M^ we see the same phenomenon but the timescale 
for such planets to become dynamically unstable is shorter, such that very few survive for 10 
Myr. We suspect that this chaotic region will be cleared out for Mg = 80 M® on timescales 
that are somewhat longer, but still short compared with the lifetime of the system. Thus, 
we do not expect any planets to exist in the 4g:ld MMR. 

4.3. Neutral Resonances — 2f:lg, 3f:lg, and 5g:2d 

4.3.1. The 2f:lg MMR 

The 2f:lg MMR is located at ~ 1.24 AU. Figure [TT] shows a lot of substructure within 
the resonance, with significant variations in FTD and stability between neighboring test 
planets. We believe these variations are caused by a combination of secular effects and 
sparse sampling. Nonetheless, we see a clear trend of higher FTD and greater stability for 
lower Mg. 

For Mg > 10 M® only a very limited sample of test planets show evidence for libration 
of 2f:lg resonant angles. Indeed, for Mg = 30 M^ and 60 M^ the only region which exhibits 
resonant libration is at = 1.24 and 1.25 AU, and Cg = 0.26-0.30. In this region libration 
of 62 = 2Xg — Xf — vof occurs but with varying amplitudes and in a chaotic fashion with 
occasional circulation. However, the median FTD value of these resonant planets is only 
0.1 {Mg = 60 M®) and 0.37 {Mg = 30 M®). A large range of parameter space exhibits 
temporary libration of resonance angles but no long-term resonance. This region is centered 
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at 1.24-1.25 with somewhat smaller eccentricities, and has small FTD values. In contrast, 
for Mg = 10 M^, several regions exhibit stable resonant libration. Resonant orbits tend to 
correlate with high FTD values in the 'V'-shaped region and tend to lie at the edges at 
= 1.24 and 1.26 AU. 

Figure [12] shows the evolution of resonant angles 6i and 62 for two simulations, both 
starting with = 1.251 AU and Cg = 0.282, but with Mg = 60 and 10 M®. For 
Mg = 60 M^, 62 librates about 0° in irregular fashion with occasional circulation, and 62 
circulates. For Mg = 10 the situation is quite different: 9i librates steadily about 75° 
with an amplitude of 30°, and 62 librates about 315° with an amplitude of ~ 90° but with 
occasional circulation!^ The contrast between the two cases is remarkable and leads us to 
the conclusion that it is very unlikely for a planet with Mg > 20 to exist in the 2f:lg 
MMR. 



4.3.2. The 3f:lg MMR 

The 3f:lg MMR lies at 1.63 AU. Figure [13] shows a clear trend between lower FTD in 
this region and larger Mg. Thus, the 3f:lg MMR is unlikely to contain a planet more massive 
than ~ 30 M^. The mean [median] values of the FTD for simulations with = 1.633 AU 
are 0.49 [0.59] for Mg = 68 M®, 0.69 [0.80] for Mg = 30 M®, and 0.97 [0.98] for Mg = 10 M®. 

None of the planets with ag = 1.633 AU in Fig. [13] (the central "column" of ag values) 
stay in resonance for long timescales. Resonant angles librate temporarily in many cases 
before switching to circulation, and sometimes back to libration in irregular fashion. Despite 
this chaotic behavior, most of these cases appear to be stable for 10 Myr, without undergoing 
close approaches with planet /. Many of the simulations with ag = 1.628 and 1.638 AU 
in Fig. [13] exhibited a period of apsidal libration between planets / and g. As for the 
resonant cases, periods of circulation and libration were often chaotically interspersed, but 
the simulations were nonetheless stable and with high FTD values. For smaller Mg, there 
exist fewer planets which exhibit temporary resonant libration, but the region of temporary 
apsidal libration is expanded. For the most part, regions of low FTD correspond to chaotic 
zones and high FTD correspond to temporary apsidal libration. 



^It is uncommon for resonant angles to librate about values other than 0° or 180° but can happen in 
some circumstances (e.g., Zhou & Sun 2003). 
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4.3.3. The 5g:2d MMR 

Figure [T^ shows the stabihty and FTD of planet g in and near the 5:2 resonance with 
planet d. The structure of the phase space is quite simple in this case and can be broken 
into four regions. The first region, represented as high-FTD areas at < 0.07, undergoes 
regular apsidal libration but is not in resonance. The second, smaller region also has high 
FTD values and is located at ag ^ 3.20 — 3.225 AU and Cg = 0.25 — 0.4. This region is 
wider for Mg = 50 than for 95 but the characteristics are the same for the two 
values of Mg-. this zone undergoes stable libration of all four resonant arguments, as well as 
apsidal libration. This region is therefore in the apsidal corotation resonance, also seen for 
large Mg in the 2g:ld MMR. The third region comprises the low-FTD region centered on the 
resonant region, at slightly smaller ag and Cg. This chaotic region is where test planets may 
be temporarily captured into resonance or apsidal libration but the evolution is chaotic and 
the resonance is short-lived. The fourth and final region includes the high-FTD areas at the 
edges of our sampled zone, at Cg > 0.1. This region does not participate in the resonance or 
apsidal libration. 

For planet g to be located in the 5g:2d MMR, it must be localized in both ag and eg. It 
must reside at ag ~ 3.21 AU with eg ~ 0.3; this resonant region is wider for lower Mg. The 
surrounding region is unlikely to host a massive planet given the low FTD values. But for 
low eg, the entire region is allowed and apsidal libration is preferred. 

4.4. The 3c:lb MMR 

Planets b and c lie very close to the 3:1 MMR (Marcy et a/. 2002; Ji et al 2003), but 
Fischer et al. (2008) note that the resonant arguments are circulating rather than librating. 
In other words, planets b and c are not in resonance. Since an additional planet g can affect 
the mean motions of other planets in the system, we calculated resonant angles of planets b 
and c for all of our stable 6-planet simulations. We find that, for our chosen configuration 
of known planets b — f, there are no cases in which planet g causes the resonant angles of 
planets b and c to librate. Thus, we conclude that the only way for planets b and c to truly 
be in a resonance is if our assumed orbital parameters for planets b — f are incorrect, which 
is certainly possible given the observational uncertainties. 
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5. Multiple Planets in the Stable Zone 

Given the width of the stable zone between planets / and d, more than one additional 
planet could exist in the region. We ran additional simulations including multiple planets 
in the stable zone. For simplicity, we chose a fixed mass of 50 for all additional planets. 
Planets were spaced such that their closest approach distances (perihelion qi vs. aphelion 
Q2) were separated by a fixed number A of mutual Hill radii Rh, where Rh = 0.5(ai + 
a2)[(Mi + M2)/3Mj^/3 (Chambers, Wetherill & Boss (1996) and subscripts 1 and 2 refer to 
adjacent planets. We ran simulations with planets spaced by A = 5 — lA.bRn in increments 
of 0.5 Rh, with five simulations for each separation with eccentricities chosen randomly to be 
less than 0.05, for a total of 100 simulations. The number of additional planets varied with 
the planet spacing, from five planets in the stable zone for A = 5 to two for A = 14.5. No 
cases with five extra planets was stable, and only one case with four extra planets survived 
for 10 Myr and the evolution of that case was chaotic. However, roughly 40% (11/28) of 
cases with three additional planets survived. Typical configurations for stable simulations 
with three planets contained planets at 1.1-1.2 AU, 1.6-1.9 AU, and 2.5-2.9 AU. The vast 
majority (43/45 = 96%) of systems with two extra planets were stable for 10 Myr. These 
contained additional planets at 1.3-1.6 AU and 2.2-3.3 AU. All stable cases had very high 
FTD values (>97%). 



6. Planets Exterior to Planet d 

Figure [13 shows the survival time of test particles beyond planet d as a function of 
their semimajor axis. As expected, there is a several AU-wide region just beyond planet d 
in which low-mass planets are unstable. In this region particles' eccentricities are quickly 
excited to values that cause them to cross the orbit of planet d, resulting in close encounters 
and ejections. Farther out, there exists a narrow contiguous region of stability from 8.6 to 9 
AU, which is roughly bounded by the 4:7 and 1:2 MMRs with planet d. This stable region 
is the only difference between our results and those of Rivera & Haghighipour (2007), who 
also mapped this outer region using test particles. The difference arises from the significant 
decrease in the best-fit eccentricity of planet d, from 0.244 to 0.063. 

A plateau of stability starts at 9.7 AU and extends continuously to 30 AU, except for 
a very narrow region of instability at the 3:1 MMR with planet d at 12.3 AU. Thus, the 
innermost planet beyond planet d is likely to be located at 10 AU or beyond, although it 
could inhabit the stable zone at 8.6-9 AU. 
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7. Conclusions 

We have mapped out the region in 55 Cancri where an additional planet g might exist. 
There is a broad region of stability between known planets / and d that could contain 
a ~Saturn-mass planet (Fig. [1]). Since observations rule out a very massive planet, our 
simulations suggest that the region could easily support two or possibly even three additional 
planets. In addition, one or more outer planets could be present in the system beyond about 
10 AU. However, such distant planets would not be detectable for many years. 

We examined eight mean motion resonances in detail (see Table 2). For two of these, 
3f:lg (i.e., the 1:3 MMR between planet / and hypothetical planet g) and 4g:ld, there was 
no stable region that exhibited regular libration of resonant arguments. Therefore, these 
resonances can not contain planets in the mass range that we explored. Given the very low 
FTD values, the 2f:lg MMR is unlikely to contain a resonant planet more massive than 
~ 20 M;^. Two other MMRs, 3g:ld and 5g:2d, may contain a stable, high-FTD resonant 
planet but the location of the MMRs is constrained to a very small region of (a^, Cg) space 
which is surrounded by a chaotic region. Finally, three MMRs, 3f:2g, 2g:ld, and 3g:2d, have 
a stabilizing influence and may contain planets near or even across the collision line with 
planet / or d. Each of these MMRs contains broad regions of stable libration of resonant 
angles, although the locations of low-FTD libration can vary with Mg. We can therefore 
only weakly constrain the presence of an additional planet in one of these resonances. 

The region between planets / and d contains many MMRs which display a wide range of 
behavior. In addition to stable and unstable resonances, the behavior of resonant arguments 
is also diverse. In some regions we would expect all resonant angles to librate regularly, but 
in others only some librate. In two instances, planet g could be in the apsidal corotation 
resonance (Michtchenko & Beauge 2003; Ferraz-Mello et a/. 2003): for large Mg in the 2g:ld 
MMR at the g — d collision line (see Fig. [5]), or in 5g:2d MMR (Fig. \\M . Moreover, we also 
see cases of "asymmetric" libration in which the equilibrium angle is neither 0° or 180° (see 
Fig. ED- Even if there are no additional planets in the f — d gap, there could be an asteroid 
belt in which this diverse and exotic dynamical behavior is on display. 

55 Cancri is a critical test of the "Packed Planetary Systems" (PPS) hypothesis, which 
asserts that any large contiguous stable region should contain a planet (BR04; RB05; Ray- 
mond et al. 2006; Barnes et al. 2008). To date, two planets have been discovered in the three 
stable zones mapped out by BR04 and RB05 (in HD 74156 and 55 Cnc). Given the width of 
the stable zone between planets / and d, PPS indicates that at least one, and possibly two 
or three, more planet (s) should exist in 55 Cancri. We look forward to further observations 
of the system that may find such planets, or perhaps show evidence of their absence. Our 
results may be used to guide observers searching for planet g and beyond. 
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Table 1. Self-Consistent Dynamical Fit of 55 Cancri (Fischer et a/. 2008) 



Planet 


M sin i (Mj) 


a (AU) 


e 


± 


zu 


Tperi (JD-2440000) 


e 


0.024 


0.038 


0.263 


0.06 


156.5 


7578.2159 


b 


0.84 


0.115 


0.016 


0.01 


164.0 


7572.0307 


c 


0.17 


0.241 


0.053 


0.052 


57.4 


7547.525 


f 


0.14 


0.785 


0.0002 


0.2 


205.6 


7488.0149 


d 


3.92 


5.9 


0.063 


0.03 


162.7 


6862.3081 




1 2 3 4 5 6 

Semimajor Axis a (AU) 



Fig. 1. — The stable zone between planets / and d. White regions represent the orbital 
elements of simulations with an additional test planet that were stable for 10 Myr. Black re- 
gions were unstable. Grey regions were stable but are unlikely to contain an additional planet 
because perturbations of the other planets' orbits were too strong (see text for discussion). 
Planets b through / are labeled. 
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Table 2. Constraints on resonant planets 



Resonance Location (AU) Comments 



2i:3g 


i.Uz-i.U4 


T") „^ j_ (2 „ ^ -| AO/1 , ,„ ,1 -1 no A 1 r\oc\ ATT ht; "L, 

Kesonant nngers at i.Uz4 and i.Uo4-i.Uo9 AU. rlign- 

FTD in fingers at e,, < 0.2. 


n.zg 


i.Zo-i.ZD 


For Mg = 30 or 60 resonance is limited to tiny region 
with very small FTD. Upper limit on resonant planet is 
~ 20 M®. 


lf:3g 


1.63 


No stable planets show resonant libration. 


4g:ld 


2.35 


No stable planets show resonant libration. 


3g:ld 


2.85-2.89 


High-FTD resonant island exists at ag — 2.86 — 2.89 
AU and Cg < 0.06. Island of apsidal hbration at ag — 
2.85 - 2.88 AU and eg = 0.15 - 0.2. 


3g:ld anti^ 


2.85-2.89 


High-FTD resonant island exists at ag = 2.86 — 2.89 AU 
and Cg < 0.01. No island of apsidal libration. 


5g:2d 


3.20 


High-FTD resonant island at ag = 3.20 — 3.225 AU and 
eg = 0.25 - 0.4. 


2g:ld 


3.7-3.8 


Resonant island at Og = 3.6 — 3.85 AU and < 0.6. 


3g:2d 


4.4-4.6 


Resonant island ag — 4.4 — 4.6 AU and = 0.1 — 0.4. 



^3:1 MMR with planet d with anti-aligned longitudes of pericenter. 
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Fig. 2. — Stability and FTD of test planets in and near the 2:3 MMR with planet / (also 
called 2/ : 3g), labeled by the test planet mass. The dashed line represents the collision line 
with planet /. Formatted as in Fig. [TJ 
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Fig. 3. — Evolution of a stable simulation in the 3f:2g MMR, with planet g starting at 1.033 
AU with Cg = 0.3. Top: Eccentricities of planets g (black) and / (grey) for a 50,000 period 
of the simulation. Bottom: Evolution of resonant argument 9i. 
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Fig. 4. — Stability and FTD of test planets in and near the 2:1 MMR with planet d (also 
called 2g : Id), labeled by the test planet mass. The dashed line is the collision line with 
planet d. Formatted as in Fig. [TJ 
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Semimajor Axis a (AU) 

Fig. 5. — The stable zone of the 2g : Id MMR, with colors that correspond to which resonant 
angles are librating. White indicates no resonant libration dark grey indicates libration of 
$2 and light grey libration of 9i, 62 and vOg — zua - this configuration is calle the apsidal 
corotation resonance (ACR). Blac areas are unstable. The dashed hne is the collision hne 
with planet d. 
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Fig. 6. — Stability and FTD of test planets in and near the 3:2 MMR with planet d (also 
called 3(7 : 2d), labeled by the test planet mass. The dashed line is the collision line with 
planet d. Formatted as in Fig. [TJ 
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Fig. 7. — Stability and FTD of test planets in and near the 3:1 MMR with planet d (also 
called 3g : Id), labeled by the test planet mass in Earth masses. Formatted as in Fig. [H 
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Fig. 8 — Evolution of two simulations for the 3g:ld MMR, both with Mg = 90 M®. Left: 
Evolution of 9z (see Eqn. 2) and eccentricities Cg and for a stable resonant planet (e^ 
shifted up by 0.05 for clarity). Right: Evolution of 6i and eg, for a chaotically-evolving 
system in the resonant region. In this case, zUg and Wd started in an anti-aligncd configuration 
and libratcd about 180° for the first ~ 1.5 Myr, while the system remained in resonance. 
This system went unstable after 3.2 Myr. 
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Fig. 9. — Stability and FTD of test planets in and near the 3:1 MMR with planet d (also 
called 3(7 : Irf), but with the longitudes of pericenter of planets g and d originally in anti- 
alignment (in Fig. [7]the apses are aligned). Again, panels are labeled by the test planet mass 
in Earth masses, and formatted as in Fig. [H 
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Fig. 10. — Stability and FTD of test planets in and near the 4:1 MMR with planet d (also 
called 4(7 : Id), labeled by the test planet mass. Formatted as in Fig. [TJ 
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Fig. 11. — Stability and FTD of test planets in and near the 1:2 MMR with planet / (also 
called 1/ : 2g), labeled by the test planet mass. Formatted as in Fig. [TJ 
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Fig. 12. — Evolution of resonant argument 6i for two simulations of the 2f:lg MMR. For the 
top panel, Mg = 60 and for the bottom panel Mg = 10 M^. 
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Fig. 13. — Stability and FTD of test planets in and near the 1:3 MMR with planet / (also 
called If : 3g), labeled by the test planet mass. Formatted as in Fig. [TJ 
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Fig. 14. — Stability and FTD of test planets in and near the 5:2 MMR with planet d (also 
called 5(7 : 2d), labeled by the test planet mass. Formatted as in Fig. [TJ 
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Fig. 15. — Survival time of test particles exterior to planet d at 5.9 AU (shown with black 
circle). Test particles extended to 30 AU; all past 15 AU were stable. 



